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ABSTRACT 

Context. Astrophysical observations must be corrected for their imperfections of instrumental origin. When Charge Coupled Devices 
(CCDs) are used, their dark signal is one such hindrance. In their pristine state, most CCD pixels are 'cool', i.e. they exhibit a low, 
quasi uniform dark current, which can be estimated and corrected for. In space, after having been hit by an energetic particle, pixels 
can turn 'hot', viz. they start delivering excessive, less predictable, dark current. The hot pixels need therefore to be flagged so that 
subsequent analysis may ignore them. 

Aims. The image data of the PICARD SODISM solar telescope (Meftah et al. 2013) require dark signal correction and hot pixel 
identification. Its E2V 42-80 CCD operates at -7.2°C and has a frame transfer architecture. Both image and memory zones thus 
accumulate dark current during, respectively, integration and readout time. These two components must be separated in order to 
estimate the dark signal for any observation. This is the main purpose of the Dark Signal Model presented in this paper. 
Methods. The dark signal time series of every pixel is processed by the 'unbalanced Haar technique' (Fryzlewicz 2007) in order 
to timestamp the instants when its dark signal is expected to change significantly. In-between those, both components are assumed 
constant, and a robust linear regression vs. integration time provides first estimates and a quality coefficient. The latter serves to assign 
definitive estimates for this pixel and for that period. 

Results. Our model is part of the SODISM Level 1 data production scheme. To check its reliability, we verify on dark frames that 
it leaves a negligible residual bias (5 e~), and generates a small RMS error (25 e~ rms). We also analyze the distribution of the image 
zone dark current. The cool pixel level is found to be 4. 1 e~ • pxL 1 • s _1 , in agreement with the predicted value. The emergence rate of 
hot pixels is investigated too. It legitimates a threshold criterion at 50 e~- pxL 1 • s~'. The growth rate is found to be on average ~500 
new hot pixels per day, i.e. 4.2% of the image zone area per year. 

Conclusions. A new method for dark signal correction of a frame transfer CCD operating at only ca. — 10°C is demonstrated. It 
allows making recommendations about the scientific usage of such CCDs in space. Independently, aspects of the method (adaptation 
of the unbalanced Haar technique, dedicated robust linear regression) have a generic interest. 

Key words. Instrumentation: detectors - Methods: data analysis - Techniques: image processing 

1 . Introduction A particular motivation to improve the removal of the dark 

signal in CCD observational data comes from the fact that it is 

When recorded for an intended astrophysical investigation or for one of the firstly applied corrections. Indeed, instrumental cor- 

any other scientific exploitation, raw observational data need to rections are not mutually commutative. They must be first esti- 

be corrected for various unwanted effects of instrumental na- mated, and then applied, in a certain logical order. The restora- 

ture. In doing so, the general goal is to reach the best possible tion steps proceed normally by rewinding upstream the flow of 

accuracy and precision in view of a faithful estimation of the information so that the remediation of detector effects comes be- 

observed physical quantity. In the present paper, we propose a fore the dealing with optical issues. As a consequence of the 

method to mitigate a couple of interrelated instrumental effects early application of the dark signal correction in the sequential 

- the dark signal and its outliers, the hot pixels - which affect cleaning of the signal, a bias at this preliminary stage propa- 

every scientific utilization of Charge Coupled Devices (CCDs) gates to posterior stages of the correction process. An under- 

unless they are cooled deeply. Our method has been designed estimation of the dark signal could e.g. convert into an over- 

for frame transfer CCDs, but some aspects of it could improve estimation of the scattered light. 

the dark signal correction of full frame CCDs. . , . • . , • ... • . c , 

b The study that is reported in this paper is part of a larger 

The dark signal is an unwanted component of the recorded effort that aims at enabling the scientific exploitation of the 



datum at each CCD picture element (pixel). It occurs due to ther- PICARD-SODISM solar images ( Meftah et al. 2013| l. Therefore, 



mally generated electric currents in the detector and is, at first the development of the method is illustrated by its application on 
order, independent on the optical quantity of interest. this special set of observations. 
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1.1. Charge Coupled Devices 

With the ever advancing CCD technology and despite the re- 
cent competition by the CMOS imaging devices, CCDs have 
constantly been, and still remain, detectors of choice for scien- 
tific applications since they were invented in 1969 at the Bell 
Telephone Laboratories (Janesick et al. 1987; Smith 2010, and 
references therein). Yet, their robust basic concept has remained 
the same: photocarriers - plus some unwanted charges generated 
by spurious sources - are collected in a potential well created 
by purposely biased electrodes and/or by ion implantation that 
generate a fixed spatially-periodic pattern of pixels in the semi- 
conductor crystal. This is the integration phase. At the end of it, 
the electric potentials of the electrodes are clocked, i.e. varied 
temporally, such that the pixel pattern slides, forcing the stored 
electrical charges to transfer along the columns, and then along 
a perpendicular linear register toward an output port where the 
signal is converted into a voltage and somehow recorded. This is 
the readout phase. See also Fig.|6]for a schematic representation 
of the CCD concept. 

The unavoidable 'imperfections' of the above physical pro- 
cess and the corrections that are required to mitigate the result- 
ing limitations have been studied by all manufacturers and scien- 
tific users of CCDs {e.g., Rodricks & Venkataraman 2005 ; Burke 
et al. 2005 ), and especially at the occasion of every space instru 
ment embarking one (e.g., Defise et al.|1997 Lo & Srour|2003 



Sirianni et al.||2004| |Schou||2004| |Penquer et al.||2009| |Gilarc ' 
etal.|2010| l. 



When CCDs are included in an instrumental setup like a 
telescope, the first correction consists in subtracting the readout 
electronic offset or bias. This first step involves the sampling of 
non-physical 'underscan' pixels and is relatively unmistakable. 
The assessment of the dark signal is the next necessary calibra- 
tion step, and it can be a delicate one as we will see in the sec- 
tions below. There are nevertheless several other performance 
issues that stem from the sensing of light and its conversion to 
a digital image by a CCD camera. The list of CCD problemat- 
ics includes the calibration of the Quantum Efficiency (QE), the 
knowledge of the video gain, the monitoring of non ideal flat- 
fields, the subtraction of image residuals (R est et al.|20 02; Crisp 
"2011) and of cosmic ray hits (CRHs) ( [Hill et al.||1997| [lpatov| 
et al.|2007 ), the taking into consideration of the worsening of the 



Charge Transfer Inefficiency (CTI) (Rhodes et al.|2010 ; Baggett 
|et al.|2012| i, etc. 



1.2. Dark current 

The work that is reported in this paper addresses the dark sig- 
nal in frame transfer CCDs, and the dark noise that is associated 
to it. Note that the dark signal is an (unsolicited) signal, and the 
dark noise is the (undesired) dark signal variability, measured by 
e.g. its RMS deviation (RMSD). Furthermore, these two distinct 
notions must not be confused either with a third one, the read 
noise (RN), which originates in the readout port of the CCD and 
in the amplification stages of the analog video chain of the cam- 
era electronics. 

The dark signal that is found in a CCD pixel comes from the 
so-called 'dark current' that is accumulated during the electronic 
integration and readout times. The dark current is due to the ther- 
mal excitation of valence electrons into the conduction band and 
to the simultaneous or subsequent collection of those electrons 
into the potential well of the pixel. It is said to be 'dark' because 
it occurs even when there is no light arriving at the CCD, e.g. 
when the shutter mechanism is closed. 



The dark signal is classically believed to grow linearly with 
the CCD integration time and to simply add arithmetically to 
the other sources of signal, namely the photoelectrons induced 
by the impinging light and the electrons produced by cosmic ray 
hits. Although some departure from linearity has been evidenced 
and modeled JWidenhorn et al.||2008] |Widenhorn eTaLl[20T0l 
|Dunlap et al.|2011||2012 1, the previous linearity statements re- 
main preponderantly valid. It will be assumed and exploited in 
the proposed method. 

There are mainly three sources of dark current: the 'surface', 



the 'depletion' and the 'diffusion' dark current fWidenhorn et al. 
2002 ). The surface dark current arises at the interface between 
the silicon and the oxide, below the electrodes. To restrain it, the 
multi-pinned phase (MPP) structure and operations are usually 
implemented. For a buried n-type channel CCD, it consists in 
biasing negatively the electrode phases and in adding a pinning 
implant under some of them. This enables signal integration hav- 
ing the Si/SiC>2 interface fully inverted. The holes accumulated 
at the interface tend to curtail the associated surface dark current 
by several orders of magnitude. Yet, some surface dark current 
can still contribute to the signal, especially when the charges are 
transfered, i.e. during the frame transfer stage if any, and during 
the readout stage. 

The depletion - or 'bulk' - dark current is produced directly 
in the depletion zone, where the electrical field of the corre- 
sponding electrode determines the potential well of the pixel. 
The 'diffusion' dark current is generated in the field-free region 
below the depletion zone. In buried channel CCDs, only the two 
latter components contribute significantly, but with different pro- 
portions, to the total dark signal. 



1.3. Hot pixels 

The dark current is not uniform. Quite the reverse, CCDs exhibit 
a number of so-called 'hot pixels' that deliver spikes of dark 
current which can be orders of magnitude larger than elsewhere 
in the frame. Spatially, they are randomly distributed and appear 
as white dots in dark images. They create an extended tail in the 
dark current distribution, and this tail determines chiefly the dark 
signal non uniformity (DSNU) (e.g., |Gilard et al.|[20CJ8| |Gilard| 
|etal.|2010p . 

A hot pixel is caused by a local discrepancy with regard to 
the perfect semiconducting crystal. It can be impurities or other 
crystallographic defects incorporated at the manufacturing stage. 
However, the advent of hot pixels in flight must be attributed 
to damages caused by energetic particles, such as protons, neu- 
trons, electrons, alpha particles, heavy ions, pions, gammas, etc. 
Those can be cosmic or solar, and many are trapped in the Van 
Allen belts ( TFeynman & Gabriel ||200Q] >. PICARD happens to 
cross the inner belt at the South Atlantic anomaly (SAA) several 
times per day. This is probably the main driver for the ignition 
of hot pixels in its CCD. 

A single event has indeed enough energy to produce either a 
transient ionizing effect that appears in a single frame and gets 
loosely labeled as a 'cosmic' ray hit (CRH), whatever its origin, 
or permanent damage that henceforth makes the pixel hot. In 
the case of inverted mode operations, granted e.g. by an MPP 
architecture, the damages are believed to be mostly displace- 
ments, typically induced by non ionizing inelastic proton col- 
lisions ( |Srour et al.||2003| |Srour & Palkol|2006] |Penquer et al] 
|2009[ ). This is why the mechanisms by which energetic protons 
degrade the CTE and produce hot pixels, have been studied thor- 
oughly for many years (Hop kinson et al.|1996||Gilard et al.| 2008 
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for example). The interaction with neutrons have received atten- 
tion too ( |Chugg et al.|2003 1 



Different studies have characterized the statistical behavior 
of hot pixels, evidencing in particular that they demonstrate the 
phenomenology of random telegraph signals (RTS) (Hopkins & 
Hopkinson 1993 Hopkinson et al. 2007 1. The capacity to anneal 



them by 'baking out' the CCD, viz. heating it above the cold op- 
erating conditions, has also been investigated (e. g.,|Defise et al. 
1997] ISirianni et aLl|2"0"04"l |Polidan et aTp004; Baggett et aT 



2012 i. In parallel, p-type channel CCDs are developed with 
some preliminary success in an attempt to surpass the optimized 
radhardness of n-type channel CCDs (Marshall et al. 2004| Gow 



et al.|2 012). Note that, albeit a nuisance, the hot pixe 



s can serve 



as a diagnostic tool to correct for the CTI (Massey et al.|2010 



1.4. Technical solutions and data processing tools for limiting 
the effects of the dark current 

For a given CCD, after design solutions have been implemented 
(e.g., Bogaart et al.|2009] ) and manufacturing care observed, the 
main and only true solution to limit the harmful effects of the 
dark signal is to cool the device. The dark current is indeed 
strongly (quasi exponentially) dependent on the temperature of 

"and Eqs. {i} & ([2} 



the silicon crystal (|Widenhorn et al. 2002 



below). But even cooled, there is still always a residual dark sig 
nal contribution, due particularly to the fact that hot pixels are 
the result of impurities or crystalline defects that reveal them- 
selves in the bulk component of the dark current, which prevails 
at low temperatures. Additionally, cooling the CCD sufficiently 
and regulating its temperature precisely require both important 
resources that are not readily available to all space instruments. 
Consequently, the CCD device cannot always be cooled as much 
or regulated as well as it should, even in major projects ( |Brown 
& Davies|2003| l. In this context, regular 'bake outs' offer com- 
plementary benefits to the cold conditions of scientific acquisi- 
tions. 

Within the range of temperatures where the CCD is normally 
operated, the data must be corrected for the dark signal if its am- 
plitude perturbs the measurement and if the dark signal correc- 
tion is indeed susceptible to enable the intended scientific ex- 
ploitation or to condition the other required instrumental correc- 
tions. To this aim, the typical (as well as minimal) approach is to 
record a dark frame just before or just after the exposed image, 
both taken with the same integration time, and to subtract the 
former from the latter. This method is straightforward but the 
mentioned subtraction increases quadratically the noise in the 
image of interest by the Poisson noise of the used instantaneous 
dark frame. It therefore impacts negatively the signal to noise 
ratio (SNR). To palliate this, a master dark frame is sometimes 
generated from multiple dark frames. However, the SNR better- 
ment offered by this master dark frame is then hampered by a 
lower duty cycle for the camera. Additionally, it is often impos- 
sible to record a dark frame, and even less so the many images 
needed for a master dark frame, using the different integration 
times of all exposed images. 

Various dedicated solutions have been proposed by different 
authors (|HiU et al.|[l997l |Widenhorn eraLl[2007l |Gilard et al. 
2010||Cai et al.|2010| for example). Although this is not the pur- 
pose of the present work, we additionally mention the approach 
of |Gomez-Rodriguez eTa l. (2009); Bu rger et aL] ( |2011| l who in- 
tend to not only correct for the dark signal but to also minimize 
the associated dark noise. 

The subject of the present paper is to report about a method 
that has been developed to estimate and update from space data 



the dark signal of the frame transfer CCD onboard the PICARD- 
SODISM solar telescope. Sect. [2] presents the relevant speci- 
ficities of the considered payload and of its space operations. 
Sect.[3]presents the principle and the details of the reconstruction 
method which leads to the desired dark signal model (DSM). 
Sect. [4] analyses and discusses the products of our DSM. The 
last section summarizes the results and concludes. 



2. Description of the context 

2.1. The PICARD mission and the SODISM instrument 

The imaging device which dark signal is modeled in the present 
paper is the flight CCD of the SODISM experiment onboard the 
PICARD space mission. PICARD was successfully launched 
on 15 June 2010 into a Sun synchronous dawn-dusk orbit, 
and commissioned in flight on 12-13 October 2010. It rep- 
resents a European asset aiming at collecting solar observa- 
tions that can serve to estimate some of the inputs to Earth cli- 
mate models (Thuillier et al. 2006). The mission scientific pay- 
load consists of the SODISM imager and of two radiometers, 
SOVAP (SOlar VAriability PICARD) and PREMOS (PREcision 
MOnitor Sensor), which carry out measurements that allow es- 
timating the Total Solar Irradiance (TSI) and the Spectral Solar 
Irradiance (SSI) from the middle ultraviolet to the red. 

The Solar Diameter Imager and Surface Mapper (SODISM) 
(Mef tah et al.|2013] > acquires continuously wide-field images of 
the photosphere and chromosphere of the Sun in five narrow pass 
bands centered at 215.0, 393.37, 535.7, 607.1, and 782.2nm. It 
contributes images that can also feed SSI reconstruction mod- 
els. Further, the scientific objectives of SODISM encompass the 
probing of the interior of the Sun via helioseismic analysis of ob- 
servations in intensity on the solar disc and at the limb (Corbard 
et al.|2008}|2013| l, and via astrometric investigations at the limb. 
The latter addresses especially the spectral dependence of the ra- 
dial limb shape, and the temporal evolution of the solar diameter 
and asphericity. 

Given the high metrological ambitions of the original 
SODISM objectives, it has been recognized that its data needed 
to be properly corrected against all relevant instrumental issues. 
Such corrections are indeed required to safeguard the morpho- 
metric and photometric potential of the data. Until an optical 
aberration became manifest, the effects of the dark signal were 
expected to count among the factors that limit the most several 
of the scientific objectives. This judgment triggered the present 
study. 

The dark signal correction is specifically intended to support 
the scientific investigations by (a) rejecting the hot pixels that 
disrupt the morphology of solar features, including sunspots and 
the radial profile of the solar limb, (b) attenuating the pattern that 
stems from the residual dark signal non uniformity (DSNU), and 
particularly the striation ensuing from the memory zone compo- 
nent of the DSNU, (c) unbiasing the subsequent estimation of 
the other sources of spurious signal, such as scattered light and 
unwanted reflections ('ghosts'), in order to minimize the photo- 
metric bias that may hamper scientific investigations. This last 
objective is quite critical because the mentioned optical effects 
must be estimated outside the solar disc where the various sig- 
nal contaminations are all unknown, weak, and of comparable 
amplitudes. 

In a prior stage, the SODISM telemetry is processed, to- 
gether with ancillary information, to produce Level (L0, or NO 
for Niveau in French) FITS files, including raw image data 
and header information. The next stage addresses the instrumen- 
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tal aspects in order to generate corrected Level 1 (LI, or Nl for 
Niveau 1 in French) products and/or LI correcting procedures, 
fed by auxiliary LI data. In this context, the reported work rep- 
resents the first step of the general LI endeavor, and our dark 
signal model is meant to supply such auxiliary LI data. 

2.2. The SODISM CCD and its camera electronics 

2.2.1. CCD procurement and characterization by the 
CO ROT program 



The SODISM CCDs are devices from the CCD 42-80 series of 
E2V They have been purchased in 2001 as part of the procure- 
ment of a batch of ten devices for the COnvection ROtation 
and planetary Transits (COROT) astrometric space mission 
(Auvergne et al. 2009 ), which has been developed a few years 
before PICARD and was also managed by the French Space 
Agency (CNES). The main rationale for selecting a device of the 
42-80 series for SODISM has been the similarity of the specifi- 
cations and the project synergy with the COROT CCD program 
at the Laboratoire a" etudes spatiales et d' instrumentation en as- 
trophysique (LESIA) jBernardi et al.||2004| [Buey et al.||2005 



|Lapeyrere et al.|2006||Giiard et al.|20d8[|Gilard et al.|2010| i 

To enhance their Quantum Efficiency (QE) the CCD 42-80 
sensors are backthinned. This technology improves the sensitiv- 
ity, especially in the middle ultraviolet (MUV), near ultraviolet 
(NUV) and blue spectral ranges. QE and pixel response non uni- 
formity (PRNU) are not further reported nor discussed here as 
those issues fall outside the scope of the present paper. 

The selection of the four COROT flight CCDs happened on 
the basis of an extensive series of tests carried out by the COROT 
CCD program (Bemardi et al. 2004; Lapeyrere et al.|2006). 



Five COROT devices that had passed the above screening 
became available to the SODISM project, who additionally pro- 
cured from E2V two more CCDs with identical specifications 
as those of the former batch (see Table [TJ. Using the outcome 
of the COROT measurements and some additional character- 
ization, the SODISM project selected for flight the CCD 42- 
80-1-985 #60 (9271-18-08), henceforth referred to as simply 
'CCD#60', or SODISM 'flight CCD'. 



Table 1. List of the CCDs available to the SODISM project. 
CCD#60 - the device selected for the SODISM flight instru- 
ment - is highlighted. 



Project origin Reference 




Name 


SODISM 


9274-03-06, 


#78 


N/A 


SODISM 


9274-18-06, 


#89 


N/A 


COROT 


9271-12-07, 


#55 


Amon 


COROT 


9271-18-08, 


#60 


Mehen 


COROT 


9235-03-07, 


#62 


Weneg 


COROT 


9235-15-06, 


#67 


Apopis 


COROT 


9272-08-07, 


#71 


Iah 



2.2.2. Format and frame transfer architecture 

The COROT/PICARD CCDs have a frame transfer architecture, 
contrarily to the regular commercial devices of the CCD 42- 
80 series. Indeed, the CCD 42 series has the flexibility for full 
frame or frame transfer variants. The commercial devices are 



full frame variants only. The reason that the full frame version 
is normally made rather than the frame transfer version (which 
could of course be operated in full frame mode) is that the addi- 
tional three bus lines required up the sides of the image area to 
drive the image and store separately make the die size slightly 
larger and reduce the buttability of the array for mosaic use. 
However, E2V has made the frame transfer variant for several 
special projects, COROT/PICARD included. The frame transfer 
architecture would allow CCD operations with no - or opened - 
shutter mechanism, and/or exposing the image zone while read- 
ing out the memory zone. Neither of these two features are used 
in SODISM. Yet, the frame transfer structure does make the dark 
signal issue worse than it would be for a full frame CCD as we 
will see in Sect. [3] 

The pixels of the CCD 42-80 series are square, have a size of 
13.5x13.5 jum 2 and no anti -blooming structures. The image zone 
(IZ) is 2048 x 2048. The IZ is where the photons of interest im- 
pinge the sensor during exposure time, when the shutter is open, 
and where the signal accumulates during the electronic integra- 
tion (see Fig. [6J. The memory zone (MZ) is 2048 x 2052. It has 
four extra lines as compared to the IZ because the serial register 
mask has four lines included with it. The MZ is masked with an 
optical aluminum shield which creates constant dark conditions 
there. This is where the signal is stored during the readout phase. 

For CCD#60, E2V indicated that its vertical CTE was bet- 
ter than 0.999999, i.e Charge Transfer Inefficiency (CTI = 
1 -CTE) < 10~ 6 at beginning of life (BOL), and that up to ~ 12 
lines of the IZ might be covered by the storeshield of the MZ. 
The LESIA tests showed that the saturation level of CCD#60 go 
from 77 10 3 e" to 1 1 1 10 3 e" from the center of the image zone to 
its edges. This was interpreted as an effect of the imperfect prop- 
agation of the electric potential to the electrodes at the middle of 
the CCD. The CCD amplifier gain being of the order of 4 p.V/e~ 
(See Sect. 2.2.4 below), the CCD output voltage is in the range 
[0;500]mV. 

The nominal SODISM cadence is one image per minute. 
This is also the maximal image cadence. Those 60 seconds in- 
clude sequentially: (a) channel selection via positioning of two 
filter wheels (except for dark frames), (b) CCD integration in the 
IZ, and associated shutter operations in case of solar exposure, 
(c) IZ->MZ frame transfer, (d) CCD readout of the MZ, and (e) 
onboard processing, including image formatting and compres- 
sion. 

2.2.3. A priori estimation of the dark current at the operating 
temperature 

The CCDs of the E2 V 42-80 series have a buried channel struc- 
ture and they function in the Advanced Inverted Mode Operation 
(AIMO). AIMO is an E2V technology which brings the benefit 
of the multi-phase pinned (MPP) structure and operations — i.e. 
negligible surface dark current - while limiting the ordinary re- 
duction in full well capacity that otherwise results from the MPP 
mode. 

The SODISM CCD is cooled so as to reach acceptable per- 
formance. For this, a radiator which is external to the structure 
of SODISM evacuates the heat from the CCD via a thermal link. 
To make the dark signal deterministic and to precisely maintain 
the pixel size in both the short and the long terms, the CCD is 
heated and regulated at -7.2 + 0.1 °C peak to peak (see Meftah 
|et al.|2013] and Figure[T]l. This is not very cold in comparison to 
other space experiments. At the CCD, the regulation wipes out 
orbital and seasonal variations of the radiator. Without regula- 
tion, the CCD would reach an unsteady temperature, wandering 
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in the range [-16; -14]°C. Temperature sensors are located be- 
low the Invar block supporting the CCDs. They provide temper- 
ature measurements with a precision of 2 mK and an accuracy of 
a few mK. 



410 b e--pxr 1 -s- 1 at +20°C, and 300 10 J e -pxr 1 -s" 1 at-7.2°C. 
If any, the hottest pixels are therefore able to produce, during a 
typical one-second integration time, a dark signal that is four 
times the saturation level of middle CCD pixels at their nomi- 
nal operating temperature of -7.2°C. One of our results below 
indicates that this apparently does not occur. 



_ 1 1 1 1 1 1 1 
r 


■ i i i i i i 


1 1 1 !. 


Lr 







Aug.201 1 



Sep. 20 12 



Fig. 1. Evolution of the temperature of the SODISM flight CCD. 
After the cooldown period in July 2010, the CCD was first regu- 
lated at -1 1°C for a month and a half. As it was anticipated that 
this temperature could not be maintained for the whole duration 
of the mission, the setup point was subsequently raised to -7.2°C 
in mid September 2010. The short term and long term stabilities 
around that point are of the order of +0.1 °C peak to peak (P- 
P). The precision of the measurement is 2 mK, viz. a hundredth 
of the amplitude of the residual variation. It can be seen that 
two bakeout periods occurred on 15-18 June 201 1 and on 13-17 
June 2012. They allowed heating the CCD up to 21°C for two 
days, and up to 24°C for three days, respectively. Brief stages 
can be noticed during their warmup and cool down phases. 



According to the manufacturer's data-sheet the dark current 
depends on temperature as per the following Arrhenius law: 



DC(T) 
DC(293 K) 



122 T J exp 



6400 \ 



122 r'exp 



2k B T 



(1) 



where DC is the dark current at temperature T [K] and DCq 
is the dark current at 293 K, both expressed in e~ ■ pxP 1 ■ s _I 
(for example). E g is the bandgap energy of Si and amounts to 
-1.1 eV; k B is the Boltzmann constant (~ 8.62 10~ 5 eV/K). Note 
that Eq. ([TJ is not in agreement with Eq. 22 of ( Widenhorn et al. 
2002| >, reproduced here: 



DC(T) = DC , m T 3 exp J-^J + DC ,de P i T 3 ' 2 exp J (2) 

LESIA measured the typical dark current of non-hot pix- 
els of CCD#60, and found O .lOe" • pxE 1 • s 1 at 233K = 
-40°C ( |Lapeyrere et al.|2006| Table 1). Following Eq. Q, the 
SODISM flight CCD should thus exhibit a 'cool' pixel dark cur- 
rent (CPDC) around 55 e~ • pxE 1 • s _I ~ 4.8 pA.cirT 2 at 293 K 
= +20°C and around 4.4 e" • pxF 1 • s -1 at -7.2°C. If we assume 
instead that the dark current is dominated by the second term of 
Eq. ([2}, non-hot pixels should exhibit a dark current of the order 
of CPDC ~ 3.6 e" • pxE 1 ■ s 1 at -7.2°C, which is 80% of our 
other estimation derived from Eq. ([JJ). 

In their CCD 42-80 datasheet, E2V also specify the hottest 
pixels to potentially deliver up to 4e~/pxl/hour at 153 K. 
Following Eq. ([TJ, those will deliver a dark current of about 



2.2.4. Serial register and readout ports 

The serial register of the E2V 42-80 CCD is split in half. This 
allows reading out the signal at the two corresponding ports. 
During nominal SODISM operations, the CCD rows are aligned 
with the polar axis of the Sun. In this configuration, the North 
pole of the Sun is in the upper half of the image which cor- 
responds to the 'right' port of the CCD. The lower half of the 
image is then readout through the 'left' port. 

The CCD 42-80 output amplifier is designed to give very 
good noise performance at low readout rates. In the datasheet, 
E2V indicates a typical read noise (RN) of 3 e~rms (with a max- 
imum of 4e~rms) at 20 kHz using correlated double sampling 
(CDS), and a typical amplifier gain of 4.5 /iV/e~ (possibly rang- 
ing from 3 to 6 //V/e~). 

For CCD#60, prior to delivery, E2V measured a RN of 3.9 
and 4. 1 e~rms for its left and right port respectively, and they 
indicated that the amplifier gains were equal to 3.96/iV/e~ for 
the left port, and to 3.85jifV/e~ for the right port, i.e. a 2.9% 
difference in favor of the left port. 

In 2003, these gains were measured by LESIA in the course 
of the COROT screening process, and they were found to be 
equal to 4.18 and 4.06yL/V/e^ at -40°C (a 3.0% difference in 
favor of the left port), with a temperature dependence measured 
to be -820ppm/K ( |Bernardi et"aT]|2004] Table 1). They were 
later reevaluated at -40°C again ( Lapeyrere et al.| |2006), and 
found to amount to 4.33 and 4.24 /uV/e~ respectively (i.e. a 2. 1 % 
difference in favor of the left port). 

Following this latest result, but using the temperature depen- 
dence found by Bernardi et al. (2004), the CCD#60 gains are ex- 
pected to be worth 4.22 and 4.14juV/e" respectively at -7.2°C, 
with a discrepancy of ~2% in favor of the left port. The exact 
value of the gain is also a function of the bias voltages and par- 
ticularly to Vod, but this extra precision is not required in the 
sequel. 



2.2.5. Camera Electronics 

The first of the three main functions of the camera electronics 
is to supply the CCD with various bias voltages. The reset drain 
voltage (Vrd) was set at 8.42 V. The substrate voltage (Vss) was 
set at 0.0 V. The output gate voltages (Vogi & ^OG2) were set at 
-6.60 V and -5.51 V, respectively. The output drain (Vod) and 
the dump drain (Vdd) voltages were set at 21.78 V and 14. 18 V, 
respectively. 

Secondly, the camera must sequence the CCD and the shutter 
mechanism with the needed clocks. It should be noted that the 
exposure defined by the mechanical shutter occurs strictly within 
the duration of the electronic integration, which always exceeds 
the former by 0.4 s for this reason. Therefore, the integration of 
dark signal lasts longer than the integration of photo-electrons 
from the solar exposure by 400 ms. The shutter exposure dura- 
tion can be programmed to take any value between 0.5 s (min- 
imum) and 16.0 s (maximum), by increments of 0.1 s, and the 
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Table 2. Parameters of the SODISM CCD camera 



Parameter 


Value (at — 7.2°C when relevant) 


CCD type 


nZW L.L.L>4z-oU series 


Plight CCD ii3mc and reference 


IVlGllGll, V^V^JJttOU 


Ullage ZU11C IDllIlaL 


Z.UH-0 A Z,UH-0 


IVlCIIlUiy ZU11C lUIIIlaL 


v 70S? 


Pixel size 


1 "K S V 1 ^ ^ urn 2 
1 j.J A 1 J.J yilll 


L, 1 1 at dUL (beginning 01 lire ) 




Saturation level 


in 111 \ rc , 
1 1 — 1 1 1 ice 


Operating temperature in flight 


i o° -L- r\ 1 
— 1 .Z ± U.1L- 


Expected dark current of non hot pixel 


~ 4 e~ ■ pxl 1 ■ s -1 


Expected dark current of hottest pixel 


300 ke" • pxl" 1 • s" 1 


Gain of the left readout port (image bottom) 


4.22/zV/e" 


Gain of the right readout port (image top) 


4.14/zV/e" 


Overall gain of the readout chain (16 bit) 


1 1 O A FiT T /_.— , A or / A HT T 

~ 1.18 ADUi6bit/e => 0.85e /ADUi 6bit 


Overall gain of the readout chain (15 bit) 


~0.59ADU 15bit /e- => 1.70e7ADU 15bit 


Read noise (RN) 


15 to 20e~rms (increasing from BOL to EOL) 


Exposure time 


0.5-16 s 


Integration time 


Exposure time + 0.4 s 


IZ-MZ frame transfer duration 


369.46 ms 


Line readout duration 


11.05 ms 


Frame readout duration 


-22 s 



integration time goes correspondingly from 0.9 s to 16.4 s. Both 
durations are given in SODISM FITS headers. 

When the integration is over, the signal that has accumulated 
in the image zone (IZ) is transfered to the memory zone (MZ) 
in 369.46 ms. This is the frame transfer and it occurs shutter 
closed, viz. under dark conditions. Note that at this point, there 
can be CTE issues occurring, as well as some evolving contri- 
bution from surface dark current which might add a column de- 
pendent pedestal. 

Thirdly, the camera electronics amplifies and converts the - 
500 mV analog signals that must be read out from the two output 
ports of the CCD (See Sect. \2.2.2\ into a stream of digital num- 
bers. The conversion cadence being 100 kHz pixel, the CCD line 
rate is 1 1 .05 ms per row, and the total readout duration amounts 
to 22 s for a whole frame that is read from the two CCD ports 
simultaneously. The digital stream is then directed to the Picard 
Gestion Charge Utile (PGCU) computing unit, where the image 
is formed, processed, and most often, compressed. 

The gain of the preamplifier stage was fixed to 9 V/V. This is 
attenuated by the impedance adaptation by a factor 0.94, which 
leads to an effective gain of 8.46 V/V. The unipolar input range of 
the analog-to-digital converter (ADC) is 0-2.8 V and it operates 
on 16 bit. This implies a conversion factor of 2 16 /2.8 analog-to- 
digital unit (ADU) per Volt, i.e. 2.34 10~ 2 ADU///V Hence, the 
overall gain of the readout chain is: 

Gie ~ 2.34 10 _2 ADU/a£V x 8.46 x 4.2juV/e" 

*0.83ADUi6bit/e-, (3) 



i.e. G, 



1.20e"/ADUi 



with a discrepancy of about 2% between the two CCD ports. 
However, most SODISM data are coded on 15 bits, the sixteenth 
and least significant bit (LSB) being simply dropped. The gain 
is then halved and becomes: 



Gis *0.42ADU 15bit /e-, 
i.e. ~ 2.40 e"/ADU 15 bit- 



(4) 



To verify the above bottom-up estimation of the gain, we 
estimate it now by means of the 'photon transfer technique' 



( |Janesick~e t al. 1987] [Downing et al.| [2006) or, more precisely, 
by using a 'dark signal transfer technique' . It states that the ran- 
dom variables of the dark signal, of its standard deviation, and 
of the read noise are related by: 

(cr s [ADU]) 2 = G[ADU/e~] x S [ADU] + (/W[ADU]) 2 (5) 

Obviously, this is valid if the noise sources are solely shot noise 
and read noise. It does not apply to the hot pixels which exhibit 
erratic fluctuations and intermittency. When RN becomes negli- 
gible with respect to the dark Poisson noise, Eq. |5| turns into: 



G _1 [e-/ADU] 



S/cr 2 s or, 



log G » log S - 2 x log cr s 



(6) 



On 31 July 2010, ten dark full frames were acquired on 
15 bits, with a 16.4 s integration time, and within few hours. 
For each pixel, after having subtracted the offsets, we com- 
pute the median and the standard deviation and plot their two- 
dimensional histogram (Figure [2]). The read noise dominated re- 
gion and the shot noise dominated region are well distinct. 

According to Eq. (|6j, log G^ 1 can therefore be estimated by 
measuring the mode of the histogram of [log (Median) - 2 log cr] 
(Figure H). Fitting a Gaussian function to this distribution gives 
a value of 1 .685+0. 1 e~/ADUi5bit for G7i. This is hardly con- 
sistent with the bottom-up estimation of 2.40e~/ADUi5bit in 
Eq. Q. The discrepancy may be explained by the approximate 
knowledge regarding e.g. the impedances used in the bottom- 
up approach. We will henceforth use the values of the gain 
measured by our in-flight analysis, which also provides an es- 
timate for the readout noise (RN), worth 9.7 ADU 15 bit, that is 
~16e~rms, in July 2010, at BOL. This is about four times 
more than the 4e~rms value expected from the CCD alone 
(Sect. |2~2~4l >. 

In order to avoid negative values that the ADC would not 
be able to convert, a voltage bias, or offset, is added to the 
signal before it is fed to the ADC. Its value is in the [840- 
850] ADU 15 bit range for the 'left' port of the CCD, and in the 
[810-820] ADU I5 bit range for its 'right' port (Figure These 
two offset values are estimated onboard for each image by cal- 
culating the average of an underscan area. The (rounded) values 
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Series of ten dark frames observed on 31 Jul 2010 with Ti— 16.4s 



. • v.'i'.cK'v' 





Aug. 2010 Dec.2010 Apr. 2011 Aug. 2011 Dec. 2011 Apr.2012 Aug. 2012 Dec. 2012 

Fig. 4. Temporal evolution of the offsets of the left port (red di- 
amonds) and right port (blue triangles) of SODISM CCD from 
July 2010 to December 2012. The offset value is larger and the 
spread is larger for the upper time series which corresponds to 
the left port. 



1 10 100 1000 

Median of the dark signal [ADU] 



10000 



Fig. 2. Dark signal transfer plot. The median and standard de- 
viation are computed for every pixel of a series of ten dark 
frames acquired on 31 July 2010 and represented here as a 2- 
dimensional histogram. A region on the left is dominated by the 
read noise (RN), which amount to about 10 ADUi5bitrms. The 
photon shot noise dominated region has a slope of 0.5 in the log- 
arithmic representation. 

Histogram of the conversion factor 

1000F ' ' — ' ' — -I 



100 



10 




0.1 1.0 10.0 

Histogram mode at 1.68546 e-/ADU 



Fig. 3. Histogram of G7i bit , the inverse of the video gain. Each 
pixel in the shot noise dominated region of Figure [2] provides 
an estimate of Gji, ir A Gaussian function is fitted and gives 
1.685+0.1 e~/ADUi5bit for the mode. 

are later transmitted in the telemetry (TM) and appended to the 
FITS header of the L0 image products. It is noticed that the time 
series of the offsets exhibit some spread and it was demonstrated 
(not shown in the present paper) that this variability is dominated 
by an orbital modulation. This coupling may be linked with the 
observed RN excess noted above. 

As the standard deviations of the same underscan area 
are also computed, this allows estimating RN. It increases 
slowly from ~ 9ADUi5bitirns * 15e~rms at BOL, up to 
~ 1 1 ADUi5bitrms m 20 e~ rms at EOL. 

The main results of Sect. l2.2l are summarized in Table|2] 



2.3. Available in-flight dark frames 

As part of the automatic ordinary science operations, also know 
as the 'SODISM routine' (Mefta h"eTal.|2013| l, it had been fore- 
seen before launch to record one dark frame per day with 1 .4 s 
integration time. This value is close to the nominal integration 
time of most channels at BOL. But it has been recognized during 
the first weeks of flight operations that dark frames with other 
integration times needed to be acquired. Indeed, for the 535 nm 
channel that is dedicated to helioseismology (535H), the expo- 
sure time is equal to 7 s. Additionally, at the occasion of special 
campaigns, the exposures in any channel can be different than 
1 s. Finally, the expected degradation of the UV channels was 
to lead to successive increases of their exposure time, which, in 
2012, has reached 16 s (the maximum value allowed by the cam- 
era electronics) for the 215 nm channel, and 6 s for the 393 nm 
channel. 

However, PICARD is a micro-mission and some time has 
been required to modify the coded routine. This is why the inte- 
gration time of the daily dark frames had to be alternated manu- 
ally (between 0.9 s and 7.4 s) in autumn 2010 (25 Sept. -2 Dec. 
2010). Later, in February and March 201 1, special linearity cam- 
paigns have permitted to vary the dark integration time at few 
occasions. But for the rest, the dark exposure time was set to 
7.4 s until 8 June 2011 when it could be again alternated man- 
ually among two values (7.4, and 16.4 s) for three months, and 
then automatically from 25 August 201 1 onwards, among three 
values (0.9, 7.4, and 16.4 s). The dark frame program has thus 
been pretty irregular during the first year of the mission when 
other aspects of the performance were at their best (Figure[5]l. 

It should be mentioned that other dark signal images have 
been acquired ~50 times per day from 5 Aug. 2010 until 27 
Mar. 2012. However, they offer information only within a ring 
shape spanning across the solar limb, and until now, those were 
not exploited by the present analysis because we are interested 
in estimating a model of the dark signal over the whole CCD 
frame. 



3. Modeling the dark signal from in-flight data 

Once modeled, the dark signal correction consists simply in a 
subtraction since we will assume linear behavior with respect to: 
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Fig. 5. Evolution of the integration time of the full CCD dark 
frames that are available to generate the desired dark signal 
model (DSM). It can be noticed that the dark frame integration 
time pattern was not regular during the first year of the mis- 
sion, until end of August 2011, after which three dark frames 
are recorded each day, with integration time equal to 0.9, 7.4 and 
16.4 s. Before the summer of 2011, most days and many weeks 
exhibit no alternation of the dark frame integration time. 



- its addition to the (wanted and spurious) components of the 
signal, and 

- integration time proportionality. 

But as a prerequisite, before performing this subtraction, we 
need to model, viz. to be able to estimate, the CCD dark signal 
for any exposure duration and at any instant of time during the 
mission lifespan. This section describes a method that delivers 
such a dark signal model (DSM) on the basis of in-flight dark 
observations scheduled with only few different integration times. 

3.1. Image zone and memory zone dark signal components 

Regarding frame transfer CCDs, the main consideration pertains 
to the generation of dark current in both the image zone (IZ) and 
the memory zone (MZ). The dark signal is indeed the addition of 
an IZ component that is proportional to the programmable inte- 
gration time, with an MZ component that is the sum of the con- 
tributions of all physical pixels of the MZ where the quantity of 
interest sojourns during readout (see Fig.|6]l. Practically, a pixel 
of coordinates (/, j), having thus accumulated signal charges in 
row i and column j of the IZ, transfers rapidly its signal, just be- 
fore the start of readout, to the pixel of same coordinates in the 
MZ. This TZ-MZ' frame transfer takes only ~370ms and is ex- 
pected to add little extra dark signal. Once in MZ, the quantity of 
interest gets shifted i times toward the serial register. During this 
slow process, it accumulates supplementary dark signal in each 
of the i physical pixels of the MZ where it dwells for 1 1.05 ms 
while one of the foregoing CCD lines is read out. 

The above scenario is formalized in Eq. (|7]) and ([8]): 



DS u (t) = O u (t)+ds u (t), 



ds{j(t) = G ■ 



(r + 5r) . T IZ (f) + T .^ T MZ (f) 



k=0 



(7) 



(8) 
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Fig. 6. Sketch of a frame transfer CCD and scenario for the dark 
signal build up. A frame transfer CCD, such as SODISM's, has 
two zones, IZ and MZ having the (quasi-) same format. During 
the CCD integration time, the signal is generated and stored in 
the Image Zone (IZ) on the right of the figure. This is when the 
shutter can be opened and the device exposed to observe the Sun. 
This is also when the IZ component of the dark signal accumu- 
lates at the rate of Tj^(f) e~ ■ pxE 1 • s _1 . The image frame is then 
transfered to the Memory Zone (MZ). At this occasion, the sig- 
nal charges that were generated at coordinates (i, j) of the IZ are 
moved to (;, j) of the MZ. During the next and last phase, the 
readout phase, the CCD lines (represented vertically here) are 
transfered one by one to the serial register. The readout phase 
lasts 22 s, which is sufficient to accumulate a significant amount 
of extra dark signal while sojourning in the i pixels (k, j) of the 
MZ. 



DS ij(t) is the random variable of the dark measurement 
recorded at date t for the pixel of coordinates (i, j). Note that 
the available dark images have always been recorded with 
15 bit resolution. 

Oij(t) is the random variable of the offset in this pixel at this 
time. Apart from the CCD port dependence, we will hence- 
forth consider that it does not depend on the pixel as we do 
not have access to this information. Hence Oij(t) = (9(f) and 
we use the amount that is computed onboard and provided in 
the L0 FITS header for the given CCD port (see Sect. [253) . 
dsij(t) is the random variable of the dark signal at t for the 
pixel (/, j). 



where, 



- G is the gain of the video chain (see Sect. 2.2.5 and Eqs. ([3) 
and ((4])). 

- T is the programmed duration of the shutter exposure, also 
known as the exposure time. 

- T + 5T — T is the duration of the CCD integration, also 
known as the integration time. As mentioned in Sect. |2.2.5| 
ST = 0.4 s. 

- T /z is the random variable of the dark current map in the 
image zone (IZ) at date t. 

- t is the duration needed to read a CCD line, r = 1 1.05 ms. 

- T MZ is the random variable of the dark current map in the 
memory zone (MZ) at date t. 

Eq. (|8]l can then be rewritten into: 

d SiJ (f) = G ■ [T' • Y0(r) + r ■ Y*f (f)] (9) 

with, 

i 

¥*f(f) = Yrffoo (io) 



k=0 
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For every pixel, there are two unknowns, T /z and *P MZ . In 
other words, we need to reconstruct the dark signal contributions 
of both the IZ and the MZ. This is a linear regression problem 
and it can be solved if there are at least two different integration 
times T for a given state of T /z and X V MZ . 

Before detailing the process of reconstruction, we notice in 
Eq. ( fT0] > that is a sum of positive contributions and must 
consequently be growing with i. 

Let us also recall that a pixel is cool (viz. non-hot) until it 
gets hit by an ionizing particle and turns hot (see Sect. |1.3[ ). The 
dsjj(t) time series are thus expected to display temporally piece- 
wise constant behaviors. They are not necessarily always grow- 
ing with time because hot pixels may cool down even if they 
appear to never come back cool. 

3.2. Generation of the dark signal model 

3.2.1. Global concept 

The goal is to estimate at any time t, T' z (f) and *P^ z (f)> the car- 
tographies of the IZ and MZ dark signal components. Given that 
Eq (|9]i depends only on unknowns that are associated to a single 
pixel (z, j), we proceed by processing the temporal dependence 
column per column to spare computer memory, and then, pixel 
per pixel. 

First, since integration time T' =7.4s is the configuration 
that has been programmed most regularly (see Fig. B), a pro- 
visional dark signal model (DSM) will be established for this 
particular value of T'. The 'DSM@7.4s' model will result from 
the fitting of piecewise constant functions to the observed time 
series because the ignition and cool down of hot pixels - whether 
in IZ or MZ - are sudden events. 

Secondly, for every pixel, we will identify the time intervals 
during which the dark signal of this pixel - and so, both of its IZ 
and MZ components - are constant according to DSM@7.4s. 

Thirdly, for each such interval and for each pixel, we per- 
form a robust regression versus integration time. Next, this en- 
ables updating the values of T' z (f) and *F^ z (f) for the successive 
periods of constancy. 

Finally, the T' z (f) and *P^ Z (0 data cubes are resampled on a 
daily grid by means of a median filter. 

The next three sections (Sect. |3.2.2| to Sect. |3.2.4| i describe 
the formation of 'DSM@7.4s'. Its outcome is then employed in 
Sect. l3~231 

3.2.2. Step 1 - Variance stabilization of the time series at 
T'=7.4s 

The unbalanced Haar transform (UHT) that is going to be needed 
in Step 3 (Sect. 3.2.4| i requires the time series to be made ap- 
proximatively homoscedastic. This is why the dsjj(f) of Eq. (|7]i 
(with t chosen so that T'(t) - 7.4 s) are first of all preprocessed 
by a generalized Box-Cox power transform ( Box & Coxj 1964 
|Sakia|1992| >: 



ds 



BoxCox i 



(dslf s (t) + a) -1 



7.4JU-1 



A-GM(ds£ s ) 



(11) 



where GM is the geometric mean. 

Indeed, it can easily be derived from Eq. |5]) that - for the 



cool pixels - ^BoxCox i ft) "will be homoscedastic if: 



which reminds of Anscombe] ( |1948|) too. For the Jul. 2010 - 
Nov. 2012 period, the smallest ds!-J s (t) observation has been 
-178 ADUisbit, which appears in agreement with the a value 
that is required to variance-stabilize the quadratic addition of 
read noise and shot noise. The above determined Box-Cox trans- 
form will of course be less effective on hot pixels, but it a poste- 
riori appears to be sufficient for the purpose of the UHT. 

The Box-Cox preprocessed time series of CCD column 
#1341 is represented in Fig.|7]-a. 

3.2.3. Step 2 - Median filtering of the time series at T'=7.4s 

Real observational data are cluttered with various outliers. In 
particular, telemetry missing blocks (TMBs) and cosmic ray hits 
(CRHs) introduce spiky features in the time series of most pix- 
els. Both of those are fortunately easy to correct for by using a 
median filter that flags the signal when it departs from the run- 
ning median by more than 5 running <x. When this occurs, the 
information is declared erroneous and replaced by the running 
median. 



3.2.4. Step 3 - Unbalanced Haar transform of the time 
series at T'=7.4s 

The sought dark signal models, and 'DSM@7.4s' in particular, 
are anticipated to be piecewise constant (viz. staircase) func- 
tions of time. Indeed, the advent of hot pixels or their partial 
cooldown, either in image zone (IZ) or memory zone (MZ), 
should create perfect steps in the DSM@7.4s time series. If in 
the IZ, the transition will be of relatively large amplitude, while 
the ignition of a hot pixel in the MZ is expected to generate a 
much smaller jump (with a ratio of t/T' w 0.14%). This discrep- 
ancy will however be attenuated by Step 1 and it remains logi- 
cal to look for an algorithm that fits staircase functions with no 
a priori on the instant of the steps, their amplitudes, nor on their 
number. Reciprocally, no change in the dark signal should oc- 
cur outside interactions with an energetic particle and annealing 
events. 

There are several such algorithms and we selected the 
multiscale 'unbalanced Haar transform' (UHT) described by 
Fryzlewicz ( 2QQ7\ because "the jumps in the basis functions do 
not necessarily occur in the middle of their support, [...] which 



avoids the restriction of jumps occuring at dyadic locations" 
(Fryzlewicz 2007), and because its computational complexity is 
of order O(ralogn). Also, [Fryzlew icz (2007) demonstrates the 



A = 0.5 and, 
a = RN 2 x G l 



Wx 1.70 = 170 ADU 



(12) 



15 bit, 



outstanding performance of his UHT scheme (see Fryzlewicz 
|2007| Figure 1, for example). 

Variance stabilization is a prerequisite to the UHT technique, 
and this motivated Steps 1 and 2. By decomposing on a Haar-like 
basis the median-filtered time series that come out from Step 2, 
the UHT provides the desired instants of time in-between which 
no hot pixel transition has taken place for the considered pixel. 

Moreover, the UHT enables denoising and reconstructing a 
piecewise constant time series that fits the input. Fig. |7]-b repre- 
sents the temporal evolution of CCD column #1341 after Step 3. 
It can be compared with Fig. |7]-a. 

For a selected pixel that ignited in the Spring of 201 1, Fig. [8] 
shows the temporal evolution of its variance-stabilized dark sig- 
nal. The UHT staircase function is overplotted. The prevalent 
goodness of the UHT fit can be appreciated in this particular 
case. 

We have reprogrammed the UHT for the Interactive Data 
Language (IDL) on the basis of the R package that Fryzlewicz] 
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(a) After Step 1 

ill 




(b) After Step 3 



Fig. 7. Upper figure (a): evolution of the dark signal in CCD col- 
umn (image row) #1341 for an integration time T equal to 7.4 s. 
The Box-Cox transformed dark signal of this CCD column are 
represented horizontally (after Step 1), and the information from 
all frames having T'=7.4s has been stacked vertically. The pres- 
ence or advent of hot pixels can be seen as bright semi-infinite 
vertical lines, telemetry holes as black horizontal segments, cos- 
mic ray hits as white dots. It can be verified that the pixels on 
the right part of the figure exhibit higher dark current since they 
dwell for a longer duration in the memory zone. We notice that 
in this particular column #1341, a hot pixel in the MZ generates 
a dark signal excess affecting the last ~10% of the CCD column, 
on the right of the figure. Lower figure (b): the same data, but 
after Step 3. 



(2007 ) makes available online A thresholding of the unbalanced 
Haar-like wavelet coefficients is needed to not reconstruct the in- 
put completely, but to preserve only the significant steps, and to 
hence denoise. We define it as per Eq. ( p"3| ), which has been de- 
vised empirically to single out the physical signal. 



\w\ X scale 225 > 4 10 4 , 



(13) 



where w is the coefficient of the considered UH wavelet, and 
scale is the length of the shortest of its two piecewise constant 
segments. Clearly, this forces short steps to display jumps of 
large amplitudes, while this permits longer periods of stability 
to result from smaller leaps. 



Variance stabilized and UHT— reconstructed time series of the dark signal 
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Fig. 8. Temporal evolution of the variance stabilized dark signal 
(T'=7.4 s) in pixel (1561,1341), over-plotted with the piecewise 
constant fit (dashed line) that results from the reconstruction by 
the unbalanced Haar technique. This represents a vertical cut in 
Fig.|7]-a and the corresponding one in Fig.|7]-b. Note that the dark 
signal appears to oscillate between two values after the pixel has 
turned hot. Such multi-level random telegraph noise (RTS) is 
expected ( |Chugg et al.|2003||Hopkinson et al.|2007 1. 



3.2.5. 



Step 4 - 
periods 



Robust regression for constant dark signal 



Having identified the instants of time when a given IZ pixel ig- 
nites, or when it changes level if it was already hot, or when 
any pixel of its corresponding MZ column discloses a signifi- 
cant jump, it becomes possible to discriminate ^ MZ and T /z , the 
MZ and IZ dark components respectively, by means of a linear 
regression on the integration time. This is the purpose of Step 4 
where we consider all dark frames that have been acquired with 
different integration times T . 

The pixel coordinates (/, j) being fixed, and t being in a time 
interval A when *P^ Z (0 and T' z (f) can be supposed constant, 
Eq. (|9| suggests that a linear regression (see e.g. |Kutner et al.| 
2004, as a possible reference) will allow separating the two. The 
regression must however be robust against the heteroscedasticity 
of the data that cannot be variance-stabilized anymore since we 
want to exploit their proportionality to T' . To this aim, we use a 
scheme that employs a weighted Li norm, as defined in Eq.(fT5]> 
below. 

We define ds^(t) = ds(t) for which t e A. Further, for each 
of the K distinct T' k , we gather the dark signals ds k A (t) defined by 
the ds&(t) for which T'(t) = 77, and we compute: 

- their median: MED k A = median,(ofs* (?)), 

- their mean absolute deviation with respect to the median: 
MAD k A = median, (\ds k A (t) - MED k A \), and 

- the maximum of all standard deviations estimated from 
Eq. @: MSD k A = max,(G ■ ds k A (t) + RN 2 ) 05 . 

As it is known that cr « 1.4826 x MeanAbsDev (Rousseeuwl 



& Leroy|2005] , we define cr k A = max[MSD*, 1.4826 • MAD*] 
in order to tentatively improve the standard deviation estimate 
for the hotter pixels. A robust linear regression (RLR) is then 
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obtained by solving iteratively: 

(T RLR , ¥f R )y = arg mini) d.s x . [ T. T I ) 

with, 



(a) CCD column #1341 Image Zone (IZ) DSM evolution 



T>0,f>() 



D 1 (ds A ,[T,¥]) = 



*6l 



(14) 



(15) 



In case K = 1, the RLR is still fed with the available ds k A (t), 
but also with ¥° SM . (the last DSM-computed value for ¥), 

Previous A r /7 

associated with a null integration time (T'=0). 



3.2.6. Step 5 - Updating and resampling on a daily grid 

We finally need to acknowledge that the intervals A are dissimi- 
lar in their capacity to inform about (Y A , *Pa)/j- For this reason, 
the (T^ 11 , ^l 11 ) that has been computed in Step 4 cannot simply 
be assigned to the period of time covered by A. The following 
updating scheme is used to adequately mix a fresh estimate com- 
puted in A by Step 4, with the previous one, already assigned to 
Previous A. 

A quality index Qa,(U) is computed for each A: 



Qi = max* T' k - min t T' k 
Qi = yfK 

Q 3 =exp(-|logD 1 (^ A ,[T A ) ¥ A ])|) 
Q A =QixQiX Q 3 



(16) 



Q\ and Q2 naturally favor the leverage provided by unalike 
and/or numerous 77 respectively. Q3 measures both the goodness 
and the reliability of the Step 4 robust linear regression (RLR). 
Note that g 3 reaches a maximum for D l (ds A , [T A ,*P A ]) = 1. 
This is the desirable behavior since the terms of the sum in 
Eq. ( p3| ) must obviously not be much larger than 1, if the fit is to 
be good. Yet, they must not be much smaller than 1 either, as this 
hints at poor statistics due to a lack of data. In case K — 1, <2a 
takes its previous value, arbitrarily divided by 50. 

Now, T and *P can be updated by taking the weighted mean 
of their new and old estimations, with their respective quality 
index as coefficients. The quality index is also updated, but using 
the geometric mean. 



•y>DSM 

A 



- (o ■ T RLR + DSM ■ T DSM 

(\cA a *^ Previous A Previa 



DSM 



Q 



DSM 



(<2a 



11/ RLR 

T A 



a) / (2a + Gp^oksa) 

0DSM . vf/DSM \ 1 1 q + /)DSM \ 

^Previous A Previous A/ ' \*^A ^-Previous A) 



DSM 

Previous A 



Qa 



(17) 

Finally, in order to provide a straightforward dark signal cor- 
rection scheme, daily IZ and MZ matrices are computed pixel 
per pixel by taking the daily median of the updating stage out- 
puts if any, or by duplicating the values of the day before oth- 
erwise. Fig. [9] shows the outcome of Step 5 for CCD column 
#1341, the same column as in Fig. [7] Fig.[T0|shows the outcome 
of Step 5 for the whole IZ and MZ as of 3 Nov. 2012. Note that 
at this stage, a binary mask of the hot pixels is additionally con- 
structed by thresholding the IZ cartographi es w ith a value (viz. 



4.2 



50 e~- pxl • s ) that is legitimated in Sect. 

The resulting daily FITS files are used by the LI process- 
ing chains of the Picard Science Mission Center (CMS-P), and 
distributed online from the SODISM homepage at LATMOS, 
which is located at http://picard.projet.latmos.ipsl. 




(b) Memory Zone (MZ) DSM evolution for the same column 




Fig. 9. Upper figure (a): evolution of the IZ dark signal model 
(DSM) for CCD column (image row) #1341. The log of the IZ 
DSM is represented horizontally, and the information from all 
frames has been stacked vertically. The advent of hot pixels in 
the IZ appears as bright semi-infinite vertical lines. Lower figure 
(b): evolution of the Memory Zone (MZ) DSM for the same col- 
umn. The representation conventions are the same as above, but 
contrarily to figure (a), the gray scale is linear. The presence or 
advent of hot pixels in the MZ is seen as steps that brighten the 
figure rightward and upward simultaneously. The dark current 
in the MZ appears as an inclined plane, which slope increases 
with time due to the gradual emergence of a distribution of hot 
pixels in the MZ. For both IZ and MZ, the first weeks (bottom 
of the images) are not well estimated because there were no data 
to feed the DSM@7.4 for those dates. 



4. Model examination 

In the present section, the products of the above described model 
are investigated. We measure in Sect. 4.1 the discrepancy be- 
tween our dark signal model (DSM) and a number of observed 
dark frames, and we verify that the bias is small. Secondly, we 
study in Sect. |4.2| the temporal evolution of the estimated IZ dark 
current. A few additional verifications are reported at the end of 
this section. 



4. 1. Checking the dark signal model (DSM) on dark frames 

The goal is here to inspect the outcome of the model in the situ- 
ation for which the solution that the DSM should ideally deliver 
is known. We select dark frames programmed with the 7.4 s in- 
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(a) Image zone dark current, 3 Nov. 2012 




(b) Memory zone dark signal, 3 Nov. 2012 




Fig. 10. Upper figure (a): DSM dark current in the image zone 
as of 3 Nov. 2012. The gray scale is logarithmic. Black corre- 
sponds to 2 ADU/s and white to 30 ADU/s. At the bottom of the 
image, one can notice an extended contamination by a spurious 
signal. Its level is in the 6-20ADU-pxL' • s _1 range (see also 
Fig. 12 1. A persistent solar image is also apparent in the IZ dark 
current. Its amplitude is about 0.2 ADUpxL 1 • s _1 . Lower figure 
(b): dark signal in the memory zone reconstructed by the DSM 
as of 3 Nov. 2012. The gray scale is linear and displays the quasi 
linear increase of the dark signal in MZ (from left to right in 
this representation). Black corresponds to 20 ADU and white to 
700 ADU. 



tegration time, i.e. the most common configuration. It is the in- 
tegration time of the images recorded in the 535 H channel that 
is dedicated to the helioseismologic investigations which are de- 
manding in term of photometric correction. This also helps get- 
ting a regular sampling, set weekly hereafter. 



(a) Evolution of the histograms of 7.4 s corrected dark frames 
1 000 I B ' ; 1 
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Fig. 11. Upper figure (a): Temporal evolution of corrected dark 
frame histograms. Dark frames having weekly regularity and 
r'=7.4s integration time have been selected. Their histogram 
is coded with a logarithmic gray scale. Most pixels display a 
corrected signal near 0, as expected. Lower figure (b): Average 
histogram of all above dark frames. Both axis have log scales. 
The number of occurrences corresponding to negative signals 
are represented by the dash line. The fitted Gaussian function 
is overplotted with a blue dotted line. Its center (the model 
bias), its standard deviation (model error), and its full width at 
half maximum (FWHM) are worth -3.1 ADU -14.7 ADU, and 
~34.6 ADU respectively. 



All along the mission, the corrected dark signal appears to 
be tightly distributed around zero, as it should (Fig. [TT[-a). The 
few lighter vertical strips are related to frames having many cos- 
mic ray hits (CRHs) taken e.g. while crossing the South Atlantic 
anomaly (SAA), and obviously, the DSM does not model the 
CRHs. By fitting a Gaussian function to the histogram (Fig.fTT} 
b), we learn that the DSM underestimates the dark signal by 
~3.1 ADU «5 e~. The model bias is thus small. 
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The Gaussian fit provides also the RMS deviation (RMSD) 
of our DSM. It is worth 14.7 ADU rms «25e"rms, i.e. the 
FWHM is 34.6 ADU -60 e", on average from Oct. 2010 till 
Oct. 2012. This is only 50% more than the read noise (RN), 
which is in the range 9-11 ADU rms for the whole period (see 
Sect. |2. 2.5 or Table|2]l. The observed RMSD accounts for possi- 
ble algorithmic inaccuracies, and certainly more predominantly, 
for phenomena that cannot be modeled, such as the RTS noise 
of hot pixels in the memory zone. 

A number of outliers depart from the Gaussian behavior and 
resulting shoulders are visible in Fig. \TT\b. By subtracting the 
above Gaussian fits from the histograms, we estimate the out- 
liers to represent ~5% of the pixels in late 2010 and ~20% in 
late 2012. These percentages cover all CRHs and hot pixels hav- 
ing non Poissonian behaviors. They cannot be estimated by a 
DSM that is based on the assumption of a sufficiently prolonged 
constancy. The error can be positive or negative. It reaches sev- 
eral thousands ADU, but this level concerns a minute fraction 
of the pixels. The solution is to flag the hot IZ pixels, so as to 
exclude them from the subsequent scientific exploitation. 

Note that columns containing hot MZ pixels cannot be alto- 
gether excluded since this would then be the case of every col- 
umn. We must consequently rely on the fact that their relatively 
unpredictable dark contributions should compensate mutually. 
However, as mentioned, this leads to a noise that probably dom- 
inates the RMSD. 



4.2. Evolution of the dark current distribution in IZ 



The normalized histogram of all DSM maps for IZ is plot- 
ted in Fig. 12 Its peak corresponds to the cool (non-hot) pixels. 
Fitting a lognormal distribution ( Baer, R. L. 2006 ), and looking 
at its mode gives 2.39 ADU-pxF • s 
dark current. Using Eq. (|4l), CPDC 



1.70e7ADU = 4.1_e 
estimated in Sect 



for CPDC. the cool pixel 
2.39 ADUpxF 1 • s 1 x 
pxF 1 ■ s~' . This is very close to the value 
which was 3.6 or 4.4 e~ • pxF 1 • s" 1 at 
— 7.2°C, depending on the physical modeling assumption. We 
remark that the persistent solar 'imprint' visible in Fig. 10 



2.2.3 



a does not show up in the histogram. It is indeed of order 
0.2 ADU pxF 1 • s" 1 (not reported in this paper) and gets merged 
in the main mode of the histogram. 

The first shoulder of the distribution, at about 6- 
20 ADU pxF 1 • s" 1 , corresponds to a pollution by an unknown 
signal near the image bottom, as already pointed out in the bot- 
tom of Fig. 10 a. 

The second shoulder is flat and extends from 30 ADU pxF 1 ■ 
s- 1 up to -2,000 ADU-pxF 1 • s" 1 * 3,400e~- pxF 1 ■ s _1 . This 
plateau represents the detectable hot pixels, which dark current 
values is quasi equiprobable over at least one order of magni- 
tude. This study allows defining a threshold at 30 ADU-pxF '-s 1 , 



viz. ~50e~ ■ pxl ■ s , for flagging the hot pixels in the image 
zone. Indeed, this value belongs to the second shoulder and not 
to the first one. Thus, the 30 ADU-pxF 1 • s _1 threshold identi- 
fies hot pixels with no risk to erroneously trigger in the bottom 
part of the image that is contaminated by signal in excess. As 
to the hottest pixels, they have a hundred times less dark current 



than specified by the manufact urer (3 ,400 e • pxl 



-l 



instead 



of 300 ke~- pxl" 1 - s"\ see Sect. 2.2.3 1. 



Hot pixels 




1000 



ADU/s 



Fig. 12. Normalized histogram of the Jul. 2010 - Nov. 2012 
maps reconstructed by the dark signal model (DSM) for the 
image zone. The histogram mode represents the cool pixels. It 
can be fitted by a lognormal function (blue dotted line), which 
gives 2.39 ADU-pxF 1 • s" 1 * 4.1 e - - pxF 1 ■ s -1 . The first small 
shoulder around 6-20 ADU-pxF 1 ■ s" 1 comprises the pixels that 
are polluted by some spurious signal in the bottom of the im- 
age. The second shoulder describes the hot pixels dark current 
distribution. It appears flat from 10-20 ADU-pxF 1 ■ s _I up to 
200 ADU-pxF 1 ■ s 1 . The hottest pixels deliver 2000 ADU-pxF 1 • 




3.5 ke • pxl 1 • s 



-l 



Oct.2010 Mar.2011 Aug. 201 1 Jan. 201 2 Jun.2012 Nov. 201 2 

Fig. 13. Temporal evolution of the fractional area covered by 
hot pixels with respect to the whole image zone. The plain line 
represents hot pixels delivering more than 250 e~pxl -s _1 , while 
the dash line corresponds to warm pixels delivering from 50 to 
250 e" ■ pxF 1 ■ s _1 . Note the two sudden drops in the otherwise 
gradual growth of the number of hotter pixels. They occur in mid 
June 201 1 and mid June 2012, when 2-3 day bakeouts happened. 



We now investigate the temporal evolution of the hot pixels. 
We split their severity in two separate ranges. The lower range 
goes from the limit of hot pixel detectivity (50 e~- pxF 1 ■ s" 1 ) up 
to 6 — 250 e" • pxF 1 ■ s _1 . The upper range contains the hotter 
pixels, delivering more dark current than 6. The value of 9 has 
been chosen to evidence the annealing threshold effect discussed 
below. The evolution of both percentages is plotted in Fig. [T3] 
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For the lower range, the rate of hot pixel ignition is slightly 
decreasing over the mission lifetime, and worth ~350 new hot 
pixels per day, viz. ~3% of the CCD image zone per year. For the 
higher range, it appears constant and worth ~150 new hot pixels 
per day, viz. ~ 1 .2% of the CCD image zone per year. This makes 
a total of ~500 new hot pixels per day, which is very comparable 
to the 1000 new hot pixels per day measured for a CCD having 
twice SODISM format by |Polidan et aT] ( |2004| >; |Baggett et IT 
( pOT2| l. Their HST/WFC3 UVIS experiment uses CCDs of the 
E2V 43 series, which format is 2048 x 4096, pixel size is 15 fim 
(instead of 13.5 /im in our case). Like SODISM's, those CCDs 
have a buried channel, an MPP implant, and are backthinned. 

We additionally remark that the 50 e~ • pxL 1 • s _I criterion 
leads to a ~ 10 % area coverage by hot pixels in late 2012, while 
we have found 20 % pixels departing from the Gaussian behav- 
ior in the previous section. This is because many of those fall 
below our threshold (see Fig.[T2"]>. But lowering it would flag too 
many pixels as being 'hot', among which many would moreover 
become false positives. 

The drops in the slope observed in mid June 2011 and mid 
June 20 1 2 in Fig. [T3] especially for the hotter pixel curve, are co- 
incident with the two SODISM bakeout operations. Indeed they 
reveal that the pixels having a dark current higher than 8 were 
instantaneously (but only partially) annealed on those occasions. 
The faster annealing of hotter pixels as compared to the anneal- 
ing of the 'warm' ones has been noted by |Polidan et al. (2004) 
and studied by Marshal l et al.| ( |2005) . An interesting possibil- 
ity inferred from the latter reference is that the hotter pixels of 
SODISM may be to some extent self-annealing constantly due 
to the 'lukecold' operating temperature of its CCD. 

Some more verifications have been made. The average dark 
current of the whole image zone is found to increase from a few 
e~-pxl~ -s _1 in mid 2010 up to SOe^-pxL'-s -1 in late 2012. This 
can be compared to the equivalent estimation that can be made 
for the memory zone. Interestingly, the integrated dark current 
in the MZ is found to be consistently 20% larger than the IZ 
counterpart. If this is real, it could mean that the shield of the 
MZ favors the generation of non ionizing collisions at the origin 
of displacement damages and hence, of hot pixels. 

Finally, it is established that the mean level of the first dark 
row of the MZ increases gradually from ~0e~/pxl in mid 2010 
up to ~40 e~/pxl in late 2012. This is possibly due to surface dark 
current seeping increasingly into the pixels during the frame 
transfer. 



5. Summary and conclusions 

Concerning the camera of the PICARD-SODISM telescope 
equipped with its Flight CCD #60, we have learned the follow- 
ing: 



1. The inverse of the camera gain is Gjg 1 « 0.85 e /ADUi6bit 
when the image is coded on 16 bits, and G^ 1 ~ 
1.70e~/ADUi5bit when it is coded on 15 bits, as is more 
commonly the case. That value has been obtained from flight 
data by means of the 'photon transfer technique' pane sick 
|et al.| 1987[ ) applied to a special sequence of dark frames 
recorded in July 2010. 

2. The overall camera read noise (RN) is measured to be 15 e~ 
rms in July 2010 and 20 e~ rms in late 2012. This is respec- 
tively four and five times larger than the intrinsic CCD read 
noise. The camera RN appears contaminated by an orbital 
variability. 



3. At the operating temperature of the Flight CCD (-7.2°C), the 
dark current of its non hot pixels is predicted to be CPDC - 
3.6 or 4.4 e~ • pxl • s _1 depending on the adopted physical 
model. 

4. The dark signal components of the CCD image zone and 
memory zone have been estimated daily by the dark signal 
model (DSM) presented in this paper. Its outcome is online 
and exploited by SODISM Level 1 products. 

5. According to our DSM, the dark current of non hot pixels 
is estimated to be CPDC - 4.1 e - pxL 1 ■ s _1 in the image 
zone, in good agreement with the expectation for a CCD at 
-7.2°C. 

6. According to our DSM, the hottest dark pixels produce 
3.5ke--pxL'- s- 1 at -7.2°C. 

7. The DSM is shown to exhibit a global bias of ~5 e~ and an 
rms deviation (RMSD) of 25 e~ rms, for images taken with a 
7.4 s integration time. 

8. A threshold can be reliably set at 50 e~ ■ pxL 1 • s _1 to flag the 
hot pixels of the image zone for their subsequent (optional) 
dismissal by science investigations. 

9. With this threshold, -1.5 % of the Flight CCD pixels are de- 
clared hot at commissioning (Oct. 2010) and ~11 % in Dec. 
2012. 

In conclusion, more general observations and results are re- 
capitulated: 



1. 



Dark signal correction is complicated when the CCD has a 
frame transfer architecture, and especially when its temper- 
ature is not cold enough and/or when the readout is not fast 
enough to essentially suppress the dark signal contribution 
of the memory zone. 

Frame transfer CCDs should therefore be employed when 
their advantages (~100% observational duty cycle, avoid- 
ance of a mechanical shutter) surpass their inconveniences. 
Temperatures below -10°C look consequently desirable for 
frame transfer CCD operations, although such 'lukecold' 
temperatures may mitigate CTE issues and the growth rate 
of hot pixels. 

It is nevertheless possible to model and correct for the dark 
signal of a frame transfer CCD by applying the method pre- 
sented in this paper. 

Aspects of our method can be useful to other applications. 
Modeling the dark signal of full frame CCDs is a natu- 
ral example. The unbalanced Haar transform of Fryz lewicz] 



2007[ ), or adapted versions thereof, appear particularly pow- 



erful to fit piecewise constant functions in various situations. 

6. Onboard estimations of e.g. the CCD offset or its read noise 
should not be rounded to integer values. 

7. It is a good idea to plan flight operations that cycle integra- 
tion time across its admissible range, automatically and reg- 
ularly, especially for the dark frames. 

8. Bakeout operations at room temperature for few days do an- 
neal a fraction of the hot pixels, and particularly the hot- 
ter ones, but loosing few days of observation can be con- 
sidered worthless with respect to their limited effectiveness. 
Reciprocally, many more bakeouts could be programmed so 
as to maintain the CCD near to its pristine state. 
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